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Abstract 

Heavy quark potential is calculated in the framework of Gauge/Gravity duality using the large D ap- 
proximation, where D is the number of dimensions transverse to the flux-tube connecting a quark and an 
anti-quark in a flat D+2 dimensional spacetime. We find that in the large D limit the leading correction to 
the ground state energy, as given by an effective Nambu-Goto string, arises not from the heavy modes but 
from the behavior of the massless modes in the vicinity of the quark and the anti-quark. This correction 
modifies the Arvis's formula for the ground state energy of a Nambu-Goto string by an overall factor which 
has a very mild dependance on the distance between the quark and the anti-quark. This analysis suggests 
that the leading corrections to an effective Nambu-Goto string description of a QCD flux-tube are in the 
form of boundary terms. 
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1 Introduction 

There is strong numerical evidence that the dynamics of a flux-tube formed between a heavy quark 
and an anti-quark is well described by the Nambu-Goto string in four space-time dimensions (see 
for e.g [1, 2, 3].) In particular the ground state energy of such a flux-tube matches very well with the 
Arvis's formula for the ground state energy of an open Nambu-Goto string with fixed end-points in 
four dimensions [4]. Yet, one expects on general grounds that there should be corrections to Arvis's 
formula, for after all, the QCD flux-tube is not a one-dimensional object, it has some "intrinsic" 
thickness which we expect to be related to the finite correlation length of the confining SU(3) 
gauge theory [5]. If one thinks of the fluctuations in the intrinsic thickness as an additional degree 
of freedom, then integrating out this degree of freedom would lead to an effective string action. 
Such an action, in general, should contain all possible terms consistent with reparameterization 
invariance. In the sprit of effective field theories, these terms can be organized into relevant, 
marginal and irrelevant terms (for a review of these ideas see, for e.g, [6, 7].) The Nambu-Goto 
action is the only relevant term, then there is a marginal term, the so called extrinsic-curvature 
term, and infinitely many irrelevant terms corresponding to higher and higher derivatives of the 
string coordinates. Further, since we are considering open strings, there could also be boundary 
terms. In fact we know on very general grounds that the resulting effective string action cannot 
be just the Nambu-Goto action [8]. 

To obtain an effective string action for QCD flux-tubes, starting from the fundamental de- 
scription in terms of quarks and gluons, requires that we integrate out all the fluctuations of the 
gauge fields whose wavelength is less than the correlation length of the confining theory. A task 
which is as difficult as to obtain an effective chiral Lagrangian for pions starting from the QCD 
action. The gauge/gravity duality [9], inspired by the AdS/CFT correspondence [10, 11, 12, 13], 
provides us with an alternate, though approximate and often heuristic, way of analyzing the dy- 
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namics of strongly coupled gauge theories. Specifically, it provides us with a geometrical way of 
exploring the consequences of one of the defining properties of SU(N) gauge theories, namely the 
existence of a finite correlation length of the gauge fields. According to the gauge/gravity duality, 
the QCD flux-tube should be thought of as a holographic projection of the fundamental string in 
five-dimensional curved space. This picture then suggests that for the confining gauge theories 
the intrinsic thickness of the flux-tube can be related to the position of a fundamental string in 
the fifth dimension [14, 15]. It was also pointed out in [16] that this relationship was consistent 
with the expected behavior of the intrinsic thickness a Nielsen- Olesen like flux-tube with its string 
tension [17]. 

The aim of the present investigation is to use the geometrical picture provided by the gauge/gravity 
duality to delineate the various sources of correction to the heavy quark potential as compared 
to the ground state energy of a Nambu-Goto string in flat four dimensional space-time. These 
corrections are induced by the fact that the fundamental QCD string lives in a five-dimensional 
curved space. We will analyze this situation using the large D approximation, where D + 2 is the 
number of dimensions of the flat spacetime where the quark lives. D is therefore the number of 
directions transverse to a QCD flux-tube. The main advantage that this approach offers is that it 
treats the massless modes of the string non-perturbatively. 

In gauge/gravity duality, the confining property of the boundary gauge theory is reflected in 
the fact that a fundamental string in five dimension can minimize its energy by staying at a fixed 
value of its fifth coordinate, Y = Y*, the coordinate which is transverse to the four dimensional flat 
Euclidian space, as shown in Fig.(l). But since the quark and the anti-quark live at the boundary, 
the end point of the fundamental string have to dip from its fixed value of Y = Y* to Y = 0, 
which is the boundary. We first calculate the large D behavior of the fundamental string in the 
approximation where the fundamental string dips from its minimal energy value to the boundary 
right at x = and x = L, where the quark and the anti-quark are placed. This is shown in Fig. (2). 
The dynamics of the open string in this approximation is very much like that of an open string 
placed at Y = Y* and stretched between two d-branes. The leading term in the large D analysis 
for this case turns out to be nothing but the Arvis's formula. As we will see in section (4), when 
one corrects for this approximation, and takes into account that the fundamental string dips not 
right at x = or x = L but there is a finite transition region marked d in Fig. (1) then the leading 
term in the large D expansion is no longer the Arvis's formula alone. Rather the Arvis's formula 
gets modified by an overall factor which can be thought of as a length dependent string tension. 

It is important to note there there should be corrections to the Arvis's formula even if the 
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Nambu-Goto action by itself provided an exact description of QCD flux-tubes [7, 18] 1 . These ex- 
pected corrections should arise from quantizing the Nambu-Goto string in non-critical dimensions. 
From the point of the large D expansion such corrections should exist since the Arvis's formula 
is only the leading term in the large D expansion and there is no a-prior reason to believe that 
the higher order corrections should vanish. Though we will not see such correction in our large 
D analysis, but it is interesting to note that these higher order corrections should themselves be 
proportional to D — 24 and therefore vanish for the critical dimension of D + 2 = 26. Thus, 
surprisingly the Arvis's formula gives the exact ground state energy of a Nambu-Goto string in 
flat space-time for D = 24 and for D — > oo. 

In an alternative approach, pioneered in [19, 20], one starts with a derivative expansion for the 
action of a fundamental string stretched between two d-branes in five dimensional confining back- 
ground, then one perturbatively integrates out the fluctuations corresponding to the oscillations 
in the fifth dimension. These world-sheet fluctuations are massive, and after integrating them out 
one obtains an effective action for the massless transverse fluctuations. In [20] it was shown that 
this leads to a correction to the Arvis's potential at the order of 1/L 4 , where L is the length of the 
fundamental string stretched between the two d-branes. This correction comes from a boundary 
term induced by integrating out the heavy modes. In our large D analysis we do not see such a 
boundary term at least to the order of 1/D, but its worth recalling that for the physical case of 
interest D = 2 and therefore even 1/D 2 correction can be important. We do see boundary terms 
but their physical origin is very different. They come not from integrating out heavy modes but 
from the behavior of the massless modes near the boundary of the five-dimensional curved space. 

The outline of the paper is as follows: In section 2, we consider a particular class of confining 
geometries and calculate the leading term in the large D expansion under the simplifying assump- 
tion that we discussed above. In section 3, within the same simplifying assumption, we calculate 
the 1/D correction arising from the fluctuation of the string in the fifth dimension which give rise 
to massive world-sheet modes. In section 4, we make allowance for the fact that the fundamental 
string reaches the boundary not suddenly at x = and x = L but smoothly over a transition region 
of finite length, d, which is a function of the distance L between the quark and the anti-quark. 
This leads to an effective string description of the QCD flux-tube in which the Nambu-Goto string 
action is supplemented by boundary terms. We summarize and state our conclusions in the final 
section. 

1 1 would like to thank Ofer Aharony for pointing this out to me. 
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2 Heavy Quark Potential from Confining Geometry in the Large D 
Limit 

The gauge invariant observable of a SU(N) gauge theory that is directly related to the heavy quark 
potential is the Wilson loop, 



Here T is a closed curve in four dimensional Euclidean space, and A is the vector potentials of 
SU(N) gauge theory [21]. According to the assumed Gauge/Gravity duality, the expectation value 
of the Wilson loop can be written as a sum over the world-sheets of an open string whose end 
points terminate on the loop T. The open string itself lives in a five dimensional curved space 
whose boundary is the the four dimensional Euclidian space where T is located [13, 22]. Formally, 



here S'fAJyjvf is the Yang-Mills action, while X(r,a) represents the world sheet of a string in a 
curved five dimensional space and S[X]ng is the Nambu-Goto action. To evaluate the Nambu- 
Goto action for a given surface in this space we need to know its metric. Ideally we would like 
a prescription which gives us the geometry of the five-dimensional space once the beta function 
(governing the running of the Yang-Mills coupling constant with energy) of the dual boundary 
quantum field theory is specified. In the absence of such a prescription we can look for simple 
modifications of AdS§ space that leads to an area law for the expectation value of the Wilson loop. 
Following [15] we will consider one such class of modification, where the metric of the curved five 
dimensional space in the coordinate system {t, x, X, Y} is given by 



In this coordinate system Y = is the four dimensional boundary 2 where the gauge theory lives 
and X denotes D = 2 directions transverse to the flux-tube lying along the x axis, while t is the 
Euclidian time. As, mentioned in the introduction, the area law for the expectation value of the 
Wilson loop is implemented by requiring that F(Y) has a minima at some fixed value of Y = Y*. 
For this class of metric we will evaluate (2.2) using a large D expansion, where D is the number 
of dimensions transverse to the QCD flux-tube. 

2 Strictly speaking for AdS§ like space, Y — is a conformal boundary [12]. F(Y) diverges there, so in what 
follows we will assume that the boundary is at Y = e and that corresponds to considering the boundary gauge-theory 
with a UV cutoff of the order 1/e [23] 




(2.1) 




(2.2) 



ds 2 = F(Y) (dt 2 + dx 2 + dX 2 + dY 2 ) . 



(2.3) 



2 Heavy Quark Potential from Confining Geometry in the Large D Limit 



6 



In the static gauge, or a physical gauge, a world-sheet of the string connecting the quark and 
the anti-quark is given by 

X = (t, x, X(t, x), Y(t, x)) = (r, a, X (a) , Y(a)) , (2.4) 

X=(X\X 2 ,...,X D ), (2.5) 

where we identify world sheet parameter r with t and a with x. A point on the world-sheet with 
coordinates (r, cr) will be represented by a. The world-sheet coordinate range from 

T T 

--<t = r<-; 0<x = a<L, (2.6) 

corresponding to a Wilson loop made up of the worldlines of a quark at x = and an antiquark 
at x = L, and we will be interested in the limit T — > oo. Surfaces appearing in (2.2) satisfy the 
boundary conditions: 

X (r, 0) = = X (r, L) ; Y (r, 0) = = Y (r, L) . (2.7) 

The minimal surface corresponding to the metric (2.3), in the limit T — > oo depends only on 
the x coordinate, 

X e =(t,x,0,Y c (x)), (2.8) 

and is qualitatively depicted in Fig.(l). Since the warp factor F(Y) has a minima at Y = Y*, 
classically the string stays at Y = Y* except at the end points where it dips to connect to the 
quark and the antiquark which are located at Y — 0. The region of transition is marked d in 
Fig.(l). For the confining geometries d is expected to grow with L but with d/L — > as L — > oo 
[24, 25], where one can define d, for example as, that distance for which 

Y*-Yc(d) _ 1A _ 3 



Y* 



1Q- S . (2.9) 



In calculating (2.2) using the large D expansion, to start with, we will neglect the region d and 
consider the fluctuation around the surface given by 

X B = {t,x,0,Y*), (2.10) 



(see Fig. (2).) In section 4 we will return to this approximation and make an estimate of the errors 
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Fig. 1: The classical string configuration in a confining geometry 




Fig. 2: Approximation to the classical string configuration used in the large D calculation 
introduced by it. 

With in this approximation the fluctuations about the minimal surface are given by 



X = (t,x,X(t,x),Y* + <f){t,x)) 



(2.11) 



with <ft(x, t) vanishing at x = and at x — L. The Nambu-Goto action for such a surface is 

"T/2 
-T/2 



Sng = T [ T/2 dt [ L dx F(Y* + <P) Vdet [ 7aft ] 

J-T/2 JO 



(2.12) 



where To is the bare string tension, and the induced world-sheet metric 7^ (apart from the warp 
factor) is 

lab = 5 ab + <9 a X ■ d b X + d a <\>d h <\>. (2.13) 
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It will be convenient to write this metric in terms of two world-sheet matrices, 

g :=g ab = S ab + d a X-d b X : (2.14) 

and 

/ := fab = d a <pd b <p. (2.15) 
Then the action for the world-sheet takes the form 

Sng = T J d 2 a F (Y* + 0) v/det^det [5 o6 + (<r7)«J- ( 2 - 16) 

One can think of T F(Y* + 0) as a position dependent string tension. Small fluctuation in the 
Y direction, 0(t, x) ^ 0, increase the string tension. These fluctuations are suppressed compared 

— * 

to the fluctuations X ^ 0, = 0, for the latter fluctuations keep the string tension fixed at its 
minimum value. This suggests the following approximate way of evaluating (2.2): Keep transverse 

— * 

fluctuations X to all orders but keep only the quadratic fluctuations in 0. In this approximation 
the Nambu-Goto action (2.12) takes the following form 

Sng 5 = \ 2 J d 2 ay/detg ab + ^ j d 2 a^detg ab Ql^f + ^M 2 2 ) (2.17) 

where M 2 is defined via, 

F (Y* + 0) = F (Y*) (l + ^M 2 2 ^) + O (0 3 ) , (2.18) 

and the fundamental length scale for the problem, l s , is given by 

1 



I 



2 T Q F(Y*). (2.19) 



In this approximation, and according to our fundamental assumption (2.2), the expectation value 
of a Wilson-loop is given by 

< W[T] > YM = J [dX] exp j-i J d 2 a^tetg~ b ^ C [g ab ] , (2.20) 

where 

C [g ab \ = J [#] exp j-i J d 2 a^detg ab QTrg^f + ^M 2 2 ^) | . (2.21) 
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Since the heavy-quark potential, V [L] , is obtained from the expectation value of a rectangular 
Wilson-Loop 

lim < W[n TxL ] > YM = exp {—TV [L]} , (2.22) 

we see that C[g a &] contains the contribution to the heavy-quark potential due to the fluctuation 
of the fundamental string in the fifth dimension. This can be heuristically interpreted as the 
contribution to the heavy-quark potential due to the fluctuation of the intrinsic thickness of the 
QCD flux-tube. 

Numerical calculations of the ground-state energy of a QCD flux-tube, in lattice gauge the- 
ory, matches very well with the ground state energy of the Nambu-Goto string in four euclidean 
dimensions as given by the Arvis's formula. This suggests that the fluctuations of the string in 
the fifth dimension should only provide a small correction to the ground-state energy given by 
the Nambu-Goto string in the four flat dimensions. Therefore, setting C[g] = constant in (2.20), 
which corresponds to describing the QCD flux-tube by a Nambu-Goto string in four Euclidean 
dimensions and ignoring the fluctuations in the intrinsic thickness, is a good approximation. 

The above observations can be made more precise by evaluating (2.20) in the large D expansion, 
where as we have noted above D is the number of transverse dimension of the D + 2 dimensional 
boundary. As we will see, in the large D expansion of (2.20) the smallness of the correction due 
to C[g] will be the consequence of the fact that there is only one massive field while there are 
D massless transverse fields, Xj. 

The large D expansion will be implemented in the standard manner [26], we start with the 
path-integral 

Z 5 = J [dX] exp j-i J d 2 oVdet0 a6 j C [g ab ] , (2.23) 

and introduce g as an independent degree of freedom in the path integral but with the constraint 
(2.14). Next we implement the constraint using the Lagrangian multiplier field N. This converts 
the integration over X fields into Gaussian integrals and one obtains 

Z 5 =Ni j [dg] [dN] exp {-S eff [g, N]} C\g], (2.24) 



where 



S e// [g,N] = i J d 2 a jVdrti+^Ok- ^} + yTr [log (-d a N ab d b )] , (2.25) 
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and Mi is a normalization constant. Further, we can integrate over the (f> field in (2.21) and write 

C[g] = A/- 2 exp [log (A)] } , (2.26) 

where 

A (g) = -d a ^tg~g' ab l d b + ^tg~M\ (2.27) 

and A/2 is another normalization constant. 

We would like to evaluate (2.23) in the limit D — > oo while keeping Dl 2 constant. In other 
words we are interested in the limit 

D ^ oo; P s = constant, (2.28) 

where I 2 is defined as 

I 2 = Dl 2 . (2-29) 
The path-integral (2.23) now can be written in the following form 

Z b =M j [dg] [dN] exp {—D S [g, N]} , (2.30) 

where S [g, N] has two parts, 

S [g, N] = ( Si [g, N] + ±-S 2 [g]) , (2.31) 



D 



with 



S 1 [g, N] = l 2 Jd 2 a{ y/deti + l -N a \5 ab - g ab ) | + ^Tr [log {-d a N ab d b )] , (2.32) 



and 



S 2 [g] = ^ 



log A(flf) . (2.33) 



Therefore in the large D limit, the contribution of the fluctuations of the field is suppressed by a 
factor of 1/D as compared to the contributions of the fluctuations of the transverse coordinates Xj. 
From the point of view of the large D expansion, the success of the four-dimensional Nambu-Goto 
string in reproducing the heavy-quark potential is due to the fact that the corrections to it are 
suppressed by a factor of 1/D. 
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The leading term in the large D expansion of (2.30) is then given by 

Z h =N h exp{-DS 1 [g,N]}, 
where the g and N are the fields that minimize the action Si alone, 

8g ' 5N 



(2.34) 



(2.35) 



Minimizing the action Si alone is noting but the large D limit of the Nambu-Goto string in a flat 
Euclidean space [26] and the required solutions are: 



g 



N 



9i 








92 


N x 








N 2 







1-A 
1-2A 

1-A 

VI - 2A 






i 



(2.36) 
(2.37) 



where A depends only on the distance between the quark and the antiquark, L, and is given by 

71 ll TlD ll 



X(L) 



24 L 2 24 L 2 



(2.38) 



A fact that will be important for us latter, when we calculate the 1/D corrections, the metric g 
has discontinuity at the boundaries, a = x = and a = x = L, which are the world-lines of the 
quark and the antiquark (see appendix of [26].) The metric at the boundary is given by 



gi> 



but interestingly N is continuous, 



N 



9bi 







1 





9b2 _ 




1-2A 



(2.39) 







VI -2A 






i 



^/l^2A 



(2.40) 



Using the above results, one readily obtains the ground state energy of a Nambu-Goto string 
in the large D limit [26], which turn out to be the Arvis's formula for the ground state energy of 
a Nambu-Goto string [4], 

DS 1 [g,N] =TV A [L], (2.41) 
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V A [L] = ^1-2X(L). (2.42) 

The large D expansion bears a close resemblance to the large N expansion of multi-component 
quantum field theories (see for e.g. [27] for a review.) It is also similar in spirit to the mean field 
approximation (for a mean field analysis of the Nambu-Goto string see [28].) Using the analogy 
with the large N expansion further, one can ask that is there an analog of the master-field for the 
SU(N) gauge theories, that is, in the large D limit can the partition function for the Nambu-Goto 
string written not in terms of sum over all surfaces, weighted by the Nambu-Goto action, but by 
a single surface. In fact there is such a fictitious mathematical surface, whose metric is given by 
the boundary metric gt, (2.39), 



J NG 



J [dX]exp|-^ j (fa^/det g al )j = J\f NG exp | - ^ j dt J dx\/det^ 
= jV JVG exp|-^ v /l-2A(L)|. 

(2.43) 

This observation will be of great use to us in evaluating the 1/D corrections in the next section. 
The metric g& should not be thought of as an induced metric, but as an independent metric defining 
the "master-surface" but which has the right boundary condition, in that when the surface meets 
the world-line of the quark and the anti-quark then g^ = 1 = g TT which is the induced metric on 
the quark and anti-quark world-lines (apart from the overall warp factor F(Y*) which is absorbed 
in our definition of l s (2.19).) 



3 The 1/D Corrections from Heavy Modes 

We will evaluate the 1/D corrections arising from the oscillation of the string along the fifth 
dimension using the observation made in the last section, namely that in the large D limit the 
partition function of the Nambu-Goto string is given by a single "master surface" with metric given 
by (2.39). 

The partition function of the fundamental string in five dimension, (2.23), can be written as 



Z 5 = Z NG < C[g] > NG = N 2 Z NG < exp <^ --Tr log ( A[; 



>NG ■ 



(3.1) 
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In the large D limit, this average in turn can be evaluated as 



TL 



1 - 2A (L) } exp <^ --Tr log A[. 



(3.2) 



where, using (2.27) and (2.39), 





d 2 2 + ^ffet^M 2 



(3.3) 



The trace can be calculated using zeta function regularization in a standard manner (see for e.g 
[29]) 



where 



and 



V„(L) 



2 



4 61 4vr 



log A 



-TV M (L) 



(3.4) 



1 9a. 

9b2 Jo 



du\og(l - e ~ 2I V^+^ , (3.5) 



1/2,, 
9b2 M - 



Using (2.39) we get 



V M (L) = --M + V h (L), 



where 



V h {L) 



4tt y/l - 2A 



dtu log ^1 



e -2L A /o; 2 +(l-2A)Af2 N j 



(3.6) 



(3.7) 



(3.f 



For LM >> 1, the integral leads to an exponential decaying term, while for LM << 1 one 
would obtain an additional 1/L correction to the heavy quark potential. Such a term is of course 
precluded both by theoretical arguments [30] and numerical results. In any case for LM << 1 our 
large D analysis fails for then the approximation of replacing F[Y C ] by F[K*] in the Nambu-Goto 
action is no longer valid, equivalently the approximation d — in Fig. (1) is no longer valid. 
Therefore the only region where the above correction is of significance is when LM ~ 1, which 
we can evaluate numerically. In Fig. (3) we show the contribution of difference between Va{L) and 
Vh(L) as compared to the Arvis's potential Va{L) for representative values of M in the units of 
inverse fermi using l s = 0.4 fermi. Therefore the main qualitative contribution of the heavy modes 
is to modify the Arvis's potential by an addition of a constant term, —\M. 
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Fig. 3: Effect of massive mode on the heavy quark potential, where AV mass = \Va — Vh\ 

4 The Boundary Effects 

Our large D calculation has been done under the approximation where the behavior of the string 
near the quark and the antiquark was ignored and effectively we have assumed that d — in 
Fig.(l). As a result of this approximation in evaluating (2.2) the fluctuations in the string world- 
sheets are incorrectly weighted near the vicinity of the quark and the anti-quark. If we relax the 
approximation d = (see Fig.(l) and Fig. (2)), then the minimal surface (2.8) is given by 



X c = [t, x, 0, Y c (x) j , 
and the fluctuations around it are now given by 

X = (t, x, X (t, x) , Y c (x) + {t, x)) . 



(4.1) 



Again, working to the quadratic order in <p (t, x), the action for these fluctuations is 

' -T/2 



S = T *J' t dt f Q dxF l Y c ( x )\ ( X + \d22 (d x Y c )^j v/deti + S h [t 



(4.2) 



where Sh[g,<p} is a quadratic functional of <i){t,x) whose explicit form is given in the Appendix 
(A.) Since the classical string configuration, Y c (x), deviates from Y = Y* only near the boundary, 
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therefore we can write the warp factor as 

F[Y C }=F* + F* V [Y C ], (4.3) 

where, to a very good approximation, r) [Y c ] has support only in the regions < x < d and 
L — d < x < L. The action then can be written as 

1 f T / 2 f L , 

S [X, 0; F (Y)} = - dt cfe^deti+ 

*>s J -T/2 JO 

\\ * f dt f dxB[Y c (x),g] v / deti+^ f ' dt [ dxB [Y c (x), g] y^tg~+ 

L s J -T/2 JO L s J -T/2 J L-d 

S h [g,<t>], (4.4) 

the world-sheet matrix g is given by (2.14) while B [Y c (x) , g] and Sh [g, 4>\ are given by (A. 4) 
and (A. 5) respectively. One can repeat the analysis of section (2), and again one finds that the 
contribution of the fluctuations <fi (t, x) along the fifth dimension, which are governed by Sh [g, 0] , 
are suppress by a factor of 1/D as compared to the mass-less fluctuations, X (t, x), which are along 
the directions transverse to the QCD flux-tube. The dynamics of these fluctuations in the large D 
limit is governed by 

"T/2 



S [X; F [Y c (x)}} = \\ dt f dx v / deti+ 

l s J -T/2 Jo 

1 f T / 2 f d , 1 f T / 2 f L , 

-/ dt dxB[Y c (x),g]y/fetfr+- dt dxB[Y c (x),g] v^diti. (4.5) 

L s J -T/2 Jo J-T/2 J L-d 

Restricting to the cases for which d/L « 1, and assuming that B[F c (a;),g] has a finite value 
at the boundary, we can approximate the last two terms by expanding the integrand about its 
boundary value to obtain 

S[X;F[Y c (x)}} = \\ ' dt f rfa;v / deti+ 

J-T/2 Jo 

- / dt (Bv^eti) +-/ dt (Bv/deti) +•••, (4.6) 

Is J-T/2 V Is J-T/2 V /x=L 



where dots represent the terms coming from derivatives of the integrand. From the point of view 
of the boundary theory, this is an effective action for a QCD flux-tube, and it is nothing but the 
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Nambu-Goto string in four flat dimensions but supplemented with boundary terms. 

Let us now come back to our assumption about B g] and inquire if it is really finite at 

the boundary 3 . From (A. 4) we have 

b[y c (x), g] = + v (*) (i + \g*2 d x rSj , (4.7) 

if the behavior of the warp factor F(Y) near the boundary is similar to that for the AdSs space 
then rj(x) near the boundary behaves as 

" (l) ~n^~P- (4 ' 8) 

where we have focused on the behavior of the fundamental string near the quark situated at the 
origin. This divergence in rj (x) is known to be related to the divergent self-energy of the quark 
and the anti-quark on which the fundamental string terminates [23] . Therefore it is plausible that 
after implementing the rernomalization procedure that removes this divergence we should be left 
with a finite boundary term of the form 

B[Y c (x), ^boundary = Z finite (g 2 2 d x Y c ) boundary = C (# 2 2 ) boundary > ( 4 - 9 ) 

where C is an unknown constant. 

If we assume that this finite boundary term can be treated as a small perturbation to the 
leading term, the Nambu-Goto term, in (4.6), then its effect can be estimated in the large D limit. 
For then the partition function in the large D limit takes the form 

Z d =j [dX] exp{-S[X;F[Y c (x)]]} 

= z ™{ e ^{-pJ_ T/ /i B ^)^-rJ_ T/ /i B ^)^)) 2 • (4 - 10) 



Evaluating it in the large D limit, using (2.43), 



T L d d 

Z d = AT d exp \ ~ pTy/frR (0; g 6 ) - —T^/g^B (L, g 6 ) \ , (4.11) 



3 I would like to thank Ofer Aharony for raising this question. 
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Fig. 4: Fluctuation leading to the formation of the caps of the flux-tube 



which leads to the following static potential 



L / d2C 



V [L] = ^ ( 1 + j — ) v^"etg fe . (4.12) 

In a confining geometry it is expected that d/L vanishes as L — > oo [24, 25], with this in mind let 
us assume that 



d I 



L VA 

where a is an unknown positive real number. Then the potential can be written as 

L ( /L\ a 1 



(4.13) 



where a and (5 are two unknown parameter which depends on the details of the confining geometry 
and A(L) is given by (2.38.) Since we haves assumed that the boundary terms are a small corrections 
to the Nambu-Goto term in (4.6) therefore we expect « 1 . Heuristically speaking, the effect 
of the boundary terms is to give the string tension of the QCD flux-tube a very mild dependence 
on the length of the flux tube. 

The need for boundary terms in a QCD flux-tube is also suggested when one considers the 
possible interaction of the QCD flux-tube with the spin of the quark and the anti-quark [31]. If 
one tries to model this interaction in an effective string description of a QCD flux-tube then one 
is forced to introduce a boundary term [32] which takes into account the intrinsic thickness of the 
QCD flux-tube. 

Finally let us consider the fluctuations depicted in Fig. (4), these are the fluctuations of the 
string in the extra dimension that have not been included in the large D calculation. In fact 
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these fluctuations cannot be described in the static gauge we have used. Though the following 
qualitative effects of these fluctuations is suggested by the idea that the holographic projection 
of a fundamental string on to the boundary is a thick flux-tube whose thickness depends on the 
distance of the fundamental string from the boundary [15]. This then implies that the shape of the 
flux-tube corresponding to the classical configuration of the fundamental string is the one shown 
in Fig. (1). The shape of the flux-tube in Fig. (2) is very different from expected shape of a QCD 
flux-tube as shown in Fig. (4), this suggests that the effect of these fluctuations in the region 
< x < d and in the region L — d < x < L is to make the flux in the flux-tube penetrate the 
region of the order of £ beyond the quark and the anti-quark. £ being the correlation length of the 
gauge field, or the inverse glue-ball mass, which one expects to be of the order of magnitude of l s . 

5 Summary and Conclusions 

The key features of the heavy quark potential as suggested by the gauge/gravity duality in the 
large D expansion can be summarized from (3.7,4.14), 

V[L] = ((1 + f(L)) + (~\m + V h [L]j , (5.1) 

where from (4.14) we see that the leading behavior of f(L) ~ 1 + (3 (^) + ■ ■ ■ ■ The contribution 
from the heavy modes is contained in V/JL] which is given by (3.8) and is shown in Fig. (3.) We 
have grouped the leading corrections in the large D expansion in the first parenthesis while the 
second parenthesis contains the 1/D corrections arising from the oscillations of the fundamental 
string in the fifth dimension. As we have noted in the introduction that there are possible 1/D 
and higher corrections arising from the massless modes in non-critical dimensions, which we have 
ignored in our analysis. There analysis and to understand the manner in which they vanish for 
the critical dimension is an important and an interesting problem to which we hope to come in 
the future. 

There are two key qualitative features of V[L]. Firstly, the correction to the Arvis's potential 
by the factor f(L). Though in our investigation we cannot fix the precise form of f(L) we can 
see clearly the reason for its existence. It arises from the simultaneous requirement that the 
fundamental string wants to lie at a fixed value Y = Y*, to minimize its energy, and the requirement 
that the fundamental string has to terminate on the boundary, Y — 0, where the quark and the 
anti-quark are placed. Thus, the origin of f(L) is in the behavior of the flux-tube in the vicinity of 
the quark and the anti-quark. Thought the present numerical results are consistent with f(L) = 0, 
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it will be very interesting to fit the future numerical calculation of the heavy quark potential to the 
above form for the heavy quark potential. The other interesting feature of (5.1) is the existence 
of the constant term in the potential, — M/4. At the first sight the existence of constant term in 
potential is of no consequence, but as has been emphasized in [33, 34], that once the renormalization 
prescription for the expectation value of the Wilson loop is appropriately fixed, the constant arises 
only from the corrections to the Nambu-Goto string 4 . It is worth recalling that there is no constant 
term in the ground state energy of a Nambu-Goto string as given by the Arvis's formula. 

The gauge/gravity duality provides us with a very useful geometrical way of understanding 
the confining property of QCD. Even though this description is largely heuristic, as we lack the 
gravitational dual of QCD that includes both confinement and the asymptotic freedom, still it 
gives us various new insight into the formation and the dynamics of the QCD flux-tube. These 
insights follow from the geometrical idea that the QCD flux-tube is a holographic projection of a 
fundamental five-dimensional string which terminates on the quark and the anti-quark that live 
on the four-dimensional boundary of a curved five-dimensional space. For example, it explains 
the success of the Nambu-Goto description of the QCD flux-tube, even though the flux-tube is 
not a one-dimensional fundamental string. From the geometrical point of view, the fluctuations of 
intrinsic thickness of the flux-tube, which distinguishes it from a fundamental string, corresponds to 
fluctuations of the fundamental string along the fifth-dimension and these fluctuations are massive 
and hence suppressed. Though the explanation again is heuristic, for the mass of these modes 
cannot be predicated and one has to appeal to dimension analysis and argue that they must be of 
the order of the glueball mass. Further the same geometric picture suggests that the Nambu-Goto 
description of the QCD flux-tube has to be supplemented by boundary terms, for in the immediate 
vicinity of the quark and the anti-quark the fluctuations of the fundamental string is qualitatively 
different from that of a Nambu-Goto string in flat four dimensional space. It is interesting to note 
that these boundary terms appear not because of integrating out the heavy mode, corresponding 
to the oscillations in fifth dimension, but because of the behavior of the massless modes near the 
boundary. Better understanding of these boundary terms therefore may also give us some clue as 
to how does gauge/gravity duality implements the asymptotic freedom. 
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A Appendix: Fluctuations Around the Minimal Surface 

We want to consider the fluctuations (4.2) around the minimal surface (2.8), their action is given 
by r2 

S=[ dt f dxF [Y c (x) + 0(x, t)} ^det (g + f c ), (A.l) 

J-T/2 JO 



' -T/2 

where g is given by (2.14), while 

(f c ) ab = d a (Y c (x) + (t, x)) d h (Y c (x) + (t, x)) . (A.2) 

If one expands the action to second order in (t, x) and use the equations of motion for Y c (x) then 
one obtains 



S [X, 0; F (Y)} = 1 / dt f dxy/fetfr + \ [ dt [ dxB [Y c (x), g] y/tet^ 

J-T/2 Jo J-T/2 Jo 

+ ^[g,0], (A.3) 



where 



D[y r (,r).g} = [ ± g ^d x Y c + r)(x) + ^ V (x)g; 2 1 d x Y c ) . (A.l) 



and 7)(x) is given by (4.3) while 
S h [g, 0] = 



dt / dx- (A [g, d x Y c ] {d t 4>f + B [g, d x Y c \ {d x <pf +) 

■T/2 JO 1 

+ / dt dx- (C [g, d x Y c \ (d t <f>) (d x <j>) + M 2 [g, d x Y c \ 2 ) . (A.5) 

J-T/2 Jo 1 



The factors A, B, C and M 2 are calculable functionals of g and d x Y but we will not require their 
explicit form. 
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